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Abstract 

This  paper  presents  issues  and  algorithms  for  the  problem  of  multiple  source  tracking  with  a 
network  of  aeroacoustic  sensors.  We  study  fusion  of  data  from  sensors  that  are  widely  separated, 
and  we  give  particular  attention  to  the  important  issues  of  limited  communication  bandwidth 
between  sensor  nodes,  effects  of  source  motion,  coherence  loss  between  signals  measured  at 
different  sensors,  signal  bandwidth,  and  noise.  We  compare  the  tracking  performance  of  vari¬ 
ous  schemes,  including  joint  (coherent)  processing  of  all  sensor  data,  as  well  as  data-reduction 
schemes  that  employ  distributed  computation  and  reduced  communication  bandwidth  with  a 
fusion  center.  Our  analysis  provides  a  quantification  of  the  potential  gain  in  source  tracking 
accuracy  that  is  achievable  with  greater  communication  bandwidth  and  joint  processing  of  sen¬ 
sor  data.  We  show  that  the  potential  gain  in  accuracy  depends  critically  on  the  scenario,  as 
determined  by  the  source  motion  parameters,  signal  coherence  between  sensors,  bandwidth  of 
the  source  signals,  and  noise  level.  For  scenarios  that  admit  increased  accuracy  with  joint  pro¬ 
cessing,  we  present  a  bandwidth-efficient  algorithm  that  involves  beamforming  at  small-aperture 
sensor  arrays  combined  with  time-delay  estimation  between  widely-spaced  sensor  arrays. 


1  INTRODUCTION 

We  are  concerned  with  tracking  moving  sources  using  a  network  of  aeroacoustic  sensors.  We 
assume  that  the  sensors  are  placed  in  an  “array  of  arrays”  configuration  containing  several  small- 
aperture  arrays  distributed  over  a  wide  area.  Each  array  contains  local  processing  capability  and  a 
communication  link  with  a  fusion  center.  A  standard  approach  for  estimating  the  source  locations 
involves  bearing  estimation  at  the  individual  arrays,  communication  of  the  bearings  to  the  fusion 
center,  and  processing  of  the  bearing  estimates  at  the  fusion  center  with  a  tracking  algorithm 
(e.g.,  see  [1,  2,  3,  4,  5]).  This  approach  is  characterized  by  low  communication  bandwidth  and  low 
complexity,  but  the  localization  accuracy  may  be  inferior  to  the  optimal  solution  in  which  the  fusion 
center  jointly  processes  all  of  the  sensor  data.  The  optimal  solution  requires  high  communication 
bandwidth  and  high  processing  complexity.  The  amount  of  improvement  in  localization  accuracy 
that  is  enabled  by  greater  communication  bandwidth  and  processing  complexity  is  dependent  on  the 
scenario,  which  we  characterize  in  terms  of  the  source  motion  parameters,  the  power  spectra  (and 
bandwidth)  of  the  signals  and  noise  in  the  sensor  data,  the  coherence  between  the  source  signals 
received  at  widely  separated  sensors,  and  the  observation  time  (amount  of  data).  We  present  a 
framework  in  this  paper  to  identify  scenarios  that  have  the  potential  for  improved  localization 


Report  Documentation  Page 


Report  Date  Report  Type 

25FEB2002  N/A 

Dates  Covered  (from..,  to) 

Title  and  Subtitle 

Contract  Number 

Issues  and  Algorithms  tor  Tracking  Multiple  Sources  with  a 
Network  of  Sensors 

Grant  Number 

Program  Element  Number 

Author(s) 

Project  Number 

Task  Number 

Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es) 

Bucknell  University  Department  of  Electrical  Engineering 
Lewisburg,  PA  17837 

Performing  Organization  Report  Number 

Sponsoring/Monitoring  Agency  Name(s)  and  Address(es) 

Department  of  the  Army,  CECOM  RDEC  Night  Vision  & 

Sponsor/Monitor’s  Acronym(s) 

Electronic  Sensors  Directorate  AMSEL-RD-NV-D  10221 

Burbeck  Road  Ft.  Belvoir,  VA  22060-5806 

Sponsor/Monitor’s  Report  Number(s) 

Distribution/ Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

Papers  from  2001  Meeting  of  the  MSS  Specialty  Group  on  Battlefield  Acoustic  and  Seismic  Sensing,  Magnetic  and 
Electric  Field  Sensors,  Volume  1:  Special  Session  held  23  Oct  2001.  See  also  ADM001434  for  whole  conference  on 
cd-rom. 

Abstract 

Subject  Terms 

Report  Classification 

unclassified 

Classification  of  this  page 

unclassified 

Classification  of  Abstract 

unclassified 

Limitation  of  Abstract 

UU 

Number  of  Pages 

19 

accuracy  relative  to  the  standard  bearings-only  tracking  method.  We  propose  an  algorithm  that  is 
bandwidth-efficient  and  nearly  optimal  that  uses  beamforming  at  small-aperture  sensor  arrays  and 
time-delay  estimation  between  widely-separated  sensors. 

The  sensor  signals  are  modeled  as  Gaussian  random  processes,  which  allows  deterministic  as 
well  as  random  propagation  effects  to  be  included.  Our  previous  work  [6,  7]  considered  a  single 
source  with  fixed  position  (no  motion).  We  extend  the  analysis  in  this  paper  to  moving  sources 
that  follow  a  parametric  motion  model. 

This  paper  is  organized  as  follows.  The  sensor  data  model  is  presented  in  Section  2  for  the 
case  of  a  non-moving  source.  Results  on  time-delay  estimation  with  partially-coherent  signals  are 
presented  in  Section  3,  which  summarizes  and  extends  the  development  in  [7].  The  sensor  data 
model  is  extended  to  moving  sources  in  Section  4.  An  algorithm  is  presented  in  Section  5,  an 
example  with  measured  aeroacoustic  data  is  included  in  Section  6,  and  concluding  remarks  are 
given  in  Section  7. 

2  DATA  MODEL  FOR  A  NON-MOVING  SOURCE 

A  model  is  formulated  in  this  section  for  the  discrete-time  signals  received  by  the  sensors  in  dis¬ 
tributed  arrays.  To  begin,  we  consider  a  single  non-nroving  source  that  is  located  at  coordinates 
(xs,ys)  in  the  (x,y)  plane.  Then  H  arrays  are  distributed  in  the  same  plane,  as  illustrated  in 
Figure  1.  Each  array  h  £  {1, . . .  ,H}  contains  Nh  sensors,  and  has  a  reference  sensor  located  at 
coordinates  (xh,Vh)-  The  location  of  sensor  n  G  {1, . . . ,  Nh}  is  at  (xh  +  Axhn,Vh  +  A yhn ),  where 
(Axfrri,  Ayhn)  is  the  relative  location  with  respect  to  the  reference  sensor.  If  c  is  the  speed  of 
propagation,  then  the  propagation  time  from  the  source  to  the  reference  sensor  on  array  h  is 

Th  =  —  =  -  \(xs  -  Xh)2  +  (ys  -  yh)2 1  /  ,  (1) 

C  C  L  J 

where  dh  is  the  distance  from  the  source  to  array  h.  We  model  the  wavefronts  over  individual  array 
apertures  as  perfectly  coherent  plane  waves.  Then  in  the  far-held  approximation,  the  propagation 
time  from  the  source  to  sensor  n  on  array  h  is  expressed  by  T}x  +  T}m •  where 

Thn  ~  —  Xs  ,  XhdAxhn  +  Vs  —A yhn  =  —  [(cos  (t>h) Axhn  +  (sin  (f>h)Ayhn]  (2) 

C  |_  (tfo  (Ifo  J  c 

is  the  propagation  time  from  the  reference  sensor  on  array  h  to  sensor  n  on  array  h,  and  (fih  is  the 
bearing  of  the  source  with  respect  to  array  h.  Note  that  while  the  far-held  approximation  (2)  is 
reasonable  over  individual  array  apertures,  the  wavefront  curvature  that  is  inherent  in  (1)  must  be 
retained  in  order  to  model  wide  separations  between  arrays. 

The  time  signal  received  at  sensor  n  on  array  h  due  to  the  source  will  be  represented  as  Sh(t—Th.— 
Thn),  where  the  vector  of  signals  s (t)  =  [si(t), . . . ,  s#(t)]T  received  at  the  H  arrays  are  modeled  as 
real-valued,  continuous-time,  zero-mean,  jointly  wide-sense  stationary,  Gaussian  random  processes 
with  —  oo  <  t  <  oo.  These  processes  are  fully  specihed  by  the  H  x  H  cross-correlation  function 
matrix 

R- sij)  =  E{s(t  +  r)  s(t)T},  (3) 

where  E  denotes  expectation,  superscript  T  denotes  transpose,  and  we  will  later  use  the  notation 
superscript  *  and  superscript  H  to  denote  complex  conjugate  and  conjugate  transpose,  respectively. 
The  ( g ,  h)  element  in  (3)  is  the  cross-correlation  function 


rs,gh{T)  =  E{sg(t  +  r)  sh(t)} 


(4) 
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Figure  1:  Geometry  of  non-nroving  source  location  and  H  distributed  sensor  arrays.  A  communi¬ 
cation  link  is  available  between  each  array  and  the  fusion  center. 

between  the  signals  received  at  arrays  g  and  h.  The  correlation  functions  (3)  and  (4)  are  equivalently 
characterized  by  their  Fourier  transforms,  which  are  the  cross-spectral  density  functions 

/OO 

rs,gh{T)exp(-ju!T)  dr  (5) 

-OO 

and  the  associated  cross-spectral  density  matrix 

Ga(u)  =  ^{R  8(t)}.  (6) 

The  diagonal  elements  GS)hh{u)  of  (6)  are  the  power  spectral  density  (PSD)  functions  of  the  signals 
s^t),  and  hence  they  describe  the  distribution  of  average  signal  power  with  frequency.  The  model 
allows  the  PSD  to  vary  from  one  array  to  another  to  reflect  propagation  differences  and  source 
aspect  angle  differences. 

The  off-diagonal  elements  of  (6),  GSi9h{ u>),  are  the  cross-spectral  density  (CSD)  functions  for 
the  signals  sg(t)  and  Sh(t)  received  at  distinct  arrays  g  /  h.  In  general,  the  CSD  functions  have 
the  form 

Gs,gh{oj)  =  Is,gh{v)  [GStgg(u)GSthh( w)]1/2,  (7) 

where  rys,gh(^)  is  the  spectral  coherence  function  for  the  signals,  which  has  the  property  0  < 
I ls,gh{w)\  <  1.  We  will  elaborate  the  interpretation  of  signal  coherence  in  Section  3,  but  co¬ 
herence  magnitude  <  1  models  random  effects  in  the  propagation  paths  from  the  source  to  arrays  g 
and  h.  Note  that  our  assumption  of  perfect  spatial  coherence  across  individual  arrays  implies  that 
the  random  propagation  effects  have  negligible  impact  on  the  intra- array  delays  r/,.n  in  (2)  and  the 
bearings  (j> i, . . .  4>h- 

The  signal  received  at  sensor  n  on  array  h  is  modeled  as  a  sum  of  the  delayed  source  signal  and 
noise, 

Zhn(t)  =  sh(t  -  rh  -  Thn)  +  whn(t),  (8) 

where  the  noise  signals  Whn{t)  are  modeled  as  real- valued,  continuous-time,  zero-mean,  jointly  wide- 
sense  stationary,  Gaussian  random  processes  that  are  mutually  uncorrelated  at  distinct  sensors,  and 
are  uncorrelated  from  the  signals.  That  is,  the  noise  correlation  properties  are 

E{Wgm{t  +  T)whn(t)}  —  r w(l~)  Sgh.Smn  (9) 

E{wgm{t  +  T)sh{t)}  =  0,  (10) 
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where  rw(r)  is  the  noise  autocorrelation  function,  and  the  noise  PSD  is  Gw( u)  =  T{rw(r)}.  We 
then  collect  the  observations  at  each  array  h  into  x  1  vectors  z ^(t)  =  [41(f), . . . ,  Zh,Nh{t)]T  for 
h  =  1 , ,H,  and  we  further  collect  the  observations  from  the  H  arrays  into  a  (N\  +  •  •  •  +  Njj)  x  1 
vector 

'  zi  (t) 

m  =  :  .  (11) 

.  z H(t)  . 

The  elements  of  Z (t)  in  (11)  are  zero-mean,  jointly  wide-sense  stationary,  Gaussian  random  pro¬ 
cesses.  We  can  express  the  CSD  matrix  of  Z (t)  in  a  convenient  form  with  the  following  definitions. 
The  array  manifold  for  array  h  at  frequency  to  is 

exp (-jujTh.i)  1  T  exp  [j^  ((cos  4>h.) Axhl  + (sin  4>h) Ay hl)] 

a h(v)  !  i  ,  (12) 

.  exp (-juTh,Nh)  J  L  exp  \j^  ((coscj)h)AxhjNh  +  (sm<j)h)AyhiNh)]  _ 

using  Thn  from  (2)  and  assuming  that  the  sensors  have  omnidirectional  response  to  sources  in  the 
plane  of  the  array.  Let  us  define  the  relative  time  delay  of  the  signal  at  arrays  g  and  h  as 

Dgh  =  Tg  —  Th.,  (13) 

where  Th  is  defined  in  (1).  Then  the  cross-spectral  density  matrix  of  Z(f)  in  (11)  has  the  form 

Gz  M  = 

(14) 

ai(w)ai(a;)HGS)ii(a;)  •••  ai(u)aH(^)H  exp(-jujD1H)GStlH(u)  ' 

•  :  +  Gw(u>)l. 

aH(u)ai(uj)H  exp(+juDiH)GS}1H(uj)*  ■■■  aH(u)aH{u:)H  Gs,Hh{u) 

The  source  CSD  functions  Gs^gh(co)  in  (14)  can  be  expressed  in  terms  of  the  signal  spectral  coherence 
7 s,gh(w)  using  (7).  Note  that  (14)  depends  on  the  source  location  parameters  (xs,ys)  through  the 
bearings  4>h  in  a/l(a;)  and  the  pairwise  time-delay  differences  Dgh . 


2.1  Cramer-Rao  Bound  (CRB) 

The  Cramer-Rao  bound  (CRB)  provides  a  lower  bound  on  the  variance  of  any  unbiased  estimator. 
The  problem  of  interest  is  estimation  of  the  source  location  parameter  vector  0  =  [xs ,  ys]T  using  T 
samples  of  the  sensor  signals  Z(0),  Z(TS), . . . ,  Z((T  —  1)  •  Ts),  where  Ts  is  the  sampling  period.  The 
total  observation  time  is  T  =  T  Ts.  Let  us  denote  the  sampling  rate  by  fs  =  1/TS  and  uis  =  27t/s. 
We  will  assume  that  the  continuous-time  random  processes  Z(t)  are  band-limited,  and  that  the 
sampling  rate  fs  is  greater  than  twice  the  bandwidth  of  the  processes.  Then  Friedlander  [8,  9]  has 
shown  that  the  Fisher  information  matrix  (FIM)  J  for  the  parameters  ©  based  on  the  samples 
Z(0),  Z(TS), . . . ,  Z((T  —  1)  •  Ts)  has  elements 


dGz(w)  i9Gz(w)  ! 

-aerGz(w>  —ae~Gzi-u) 


i,j  =  1,2, 


where  “tr”  denotes  the  trace  of  the  matrix.  The  CRB  matrix  C  =  J^1  then  has  the  property  that 
the  covariance  matrix  of  any  unbiased  estimator  0  satisfies  Cov(0)  —  C  >  0,  where  >  0  means  that 
Cov(0)  —  C  is  positive  semidefinite  [10].  Equation  (15)  provides  a  convenient  way  to  compute  the 
FIM  for  the  distributed  sensor  array  model  as  a  function  of  the  signal  coherence  between  distributed 
arrays,  the  signal  and  noise  bandwidth  and  power  spectra,  and  the  sensor  placement  geometry.  The 
CRB  is  evaluated  for  various  scenarios  in  [6,  7]. 


3  TIME-DELAY  ESTIMATION  (TDE) 

Let  us  parameterize  the  model  in  (14)  by  the  bearings  (ph  and  the  time-delay  differences  Dy Then 
we  must  address  the  issue  of  time-delay  estimation  with  signals  that  are  partially  coherent  when 
\ls,gh\  <  1-  We  consider  this  problem  first  for  the  case  of  H  =  2  sensors,  as  illustrated  in  Figure  2a 
with  the  differential  time  delay  defined  as  D  =  D21.  It  follows  from  (14)  that  the  CSD  matrix  of 
the  sensor  data  in  Figure  2a  is 


CSD 


Zl(t) 

z2  (t) 


(16) 
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e-J'“D7s,12(w)*  [G'Siii(cu)GSi22(^)]1'/2 


e+j^Dj Sjl2(a;)  [G'S)ii(w)GSi22(w)]1//" 

Gs,  22  (w) 


The  signal  coherence  function  7^12  (w)  describes  the  degree  of  correlation  that  remains  in  the 
signal  emitted  by  the  source  at  each  frequency  ui  after  propagating  to  sensors  1  and  2.  Since  the 
sensor  signals  are  modeled  as  Gaussian  random  processes,  the  coherence  loss  can  be  equivalently 
represented  in  terms  of  a  “coherent”  signal  component  s(t)  and  additional  additive  noise  processes 
n\ (t) ,  ri2  (t) ,  as  depicted  in  Figure  2b.1  The  equivalence  is 


Z\(t)  =  Si(t)  +  Wi(t)  =  (hi*s)(t)  +  ni(t)  +  Wi(t)  ,  s 

Z2(t)  =  S2(t  -  D)  +  W2(t)  =  (h2*s)(t  —  D)  +  n2(t)  +  u;2(t) 

where  *  denotes  convolution.  The  parameters  in  the  CSD  (16)  are  related  to  the  filters  Hi( u)  = 

lF{hi(t)},  i  =  1,2,  the  PSDs  Gi(u),  i  =  1,2  of  the  noise  processes2  rii(t),  i  =  1,2  and  the  PSD 
Gc( u)  of  s(t)  as  follows: 


Hi(u)  =  Gs,  nM1/2  (18) 

iLM  =  ^sMf*Gs,22^)1/2  (19) 

|7s,12H| 

PSD  [s(t)]  =  Gc(u)  =  |7s,i2M|  (20) 

PSD[m(t)]  =  Gi(cu)  =  GS)iiH[1-|7M2M|]  (21) 

PSD[n2(t)]  =  G2(u)  =  Ga)22H[l-|7M2M|].  (22) 


Note  that  we  can  define  a  “coherent”  signal-to-noise  (SNR)  ratio  at  each  sensor,  based  on  the 
coherent  signal  component  in  Figure  2b: 


SNRCii(cu) 


K12MI 


1  -  |7s,12M|  + 
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|7s,12M| 

1  -  |7s,12M|’ 


i  =  1,2. 


(23) 


The  coherent  SNR  is  severely  limited  when  there  is  coherence  loss  (|7s,i2(w)|  <  1)  between  the 
sensors,  even  if  the  source  power  is  very  large  (Gs,u(u)  — >  00). 

lrThe  equivalent  model  is  developed  in  [7],  and  it  is  also  discussed  in  the  context  of  ultrasound  image  speckle  in 
[12]- 

2The  processes  ni(t),ri2(t)  are  independent,  zero  mean,  stationary  Gaussian  random  processes  that  are  indepen¬ 
dent  from  the  signals  and  noise  wi(t) ,  W2(t) ■ 


Source, 


D 


Sensor  1 


Zjft)  =  Sjft)  +  Wift) 

„  v 

Sensor  2 

z2(t)  =  s2(t  -  D)  +  w2(t)*'''”  Additive  noise  at  sensor 


(a) 


Filters  that 
model  source 
aspect  angle  & 
deterministic 
propagation 
effects 


"Excess"  additive 
noise  causing 
partial  coherence 
(due  to  random 
propagation 
effects ) 


Additive 

sensor 

noise 


(b) 


Figure  2:  (a)  Time-delay  estimation  (TDE)  problem  for  a  non-moving  source,  (b)  Representation 
of  partial  coherence  between  si(f),  s2(f)  as  excess  additive  noise  ni(f),n2(f). 


3.1  TDE  Performance  Bounds 


In  this  section,  we  summarize  and  further  study  performance  bounds  on  time-delay  estimation 
(TDE)  with  partially  coherent  signals  that  were  originally  presented  in  [6,  7].  First,  the  CRB  on 
estimating  the  time-delay  D  in  (16)  is,  using  (15), 
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where  T  is  the  total  observation  time  of  the  sensor  data  and  SNRCij(u;),  the  coherent  SNR,  is 
defined  in  (23).  Let  us  consider  the  case  in  which  the  signal  PSDs,  noise  PSDs,  and  coherence  are 
flat  (constant)  over  a  frequency  band  from  f\  to  /2  Hz,  and  let  us  define  SNR,;  =  Gs^n/Gw,  i  =  1,2. 
Then  the  CRB  in  (24)  reduces  to3 
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3The  CRB  in  (26)  can  also  be  derived  from  results  in  [11]  based  on  the  equivalent  model  in  (18)- (22). 


(27) 


where  (27)  is  the  high-SNR  limit.  Note  that  partially  coherent  signals  |7s,i2|  <  1  limit  the  TDE 
accuracy  in  a  way  that  cannot  be  improved  by  increased  source  power.  Improved  TDE  accuracy 
is  obtained  with  partially  coherent  signals  by  increased  observation  time  T  or  increased  signal 
bandwidth  [/i,/2].  The  source  signal  bandwidth  is  not  controllable  in  passive  source  localization 
applications,  so  increased  observation  time  is  the  only  means  for  improving  the  accuracy  of  TDE 
with  partially  coherent  signals.  Source  motion  becomes  more  important  during  long  observation 
times,  and  we  add  motion  into  the  model  in  Section  4. 

It  is  well-known  that  the  CRB  on  TDE  is  achievable  only  when  the  coherent  SNR  at  the  two 
sensors  exceeds  a  threshold  [13].  For  the  case  of  TDE  with  partially  coherent  signals,  we  showed  in 
[7]  that  a  similar  threshold  phenomenon  occurs  with  respect  to  coherence.  That  is,  the  coherence 
must  exceed  a  threshold  in  order  to  achieve  the  CRB  (24)  on  TDE.  We  can  state  the  formula 
for  threshold  coherence  for  the  following  scenario.  The  signal  and  noise  spectra  are  flat  over  a 
bandwidth  of  Aw  rad/sec  centered  at  rad/sec,  and  the  observation  time  is  T  seconds.  Further, 
assume  that  the  signal  PSDs  are  identical  at  each  sensor,  and  define  the  following  constants  for 
notational  simplicity: 


G's,ii(cwo)  —  22(^0)  —  Gs,  Gw(u  0)  —  Gw,  73,12(^0)  —  Is- 


(28) 


Then  combining  the  equivalent  model  (18)-(22)  with  previously  developed  Ziv-Zakai  bounds  on 
TDE  [13],  we  can  show  that  the  following  SNR-like  expression  characterizes  the  performance  of 
time  delay  estimation  with  partially  coherent  signals: 


SNR(7s) 
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-1 
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(29) 


The  threshold  SNR  for  CRB  attainability  [13]  is  a  function  of  the  tinre-bandwidth  product  ^ 


and  the  fractional  bandwidth 
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where  (p{y )  =  l/\/27r  f°°  exp(— f2/2)  dt.  It  follows  that  the  threshold  coherence  value  is 
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(31) 


For  a  specific  narrowband  time  delay  estimation  scenario,  the  threshold  SNR  for  CRB  attainability 
is  given  by  (30),  and  (31)  provides  a  corresponding  threshold  coherence  for  CRB  attainability.  Since 
1 7s  1 2  <  1,  (31)  is  useful  only  if  Gs/Gw  >  SNRthresh. 

Figure  3  contains  plots  of  the  threshold  coherence  in  (31)  as  a  function  of  the  tinre-bandwidth 
product  (  A2rT ) ,  SNR  and  fractional  bandwidth  (-77)-  Note  that  ^  =  10  dB  is  nearly  equiv¬ 
alent  to  §r-  — »  00.  The  variation  of  threshold  coherence  with  fractional  bandwidth  is  illustrated  in 
Figure  3d.  We  note  that  very  large  tinre-bandwidth  product  is  required  to  overcome  coherence  loss 
when  the  fractional  bandwidth  is  small  at  0.1.  For  a  fixed  threshold  coherence  value,  such  as  0.7, 
each  doubling  of  the  fractional  bandwidth  reduces  the  required  tinre-bandwidth  product  by  about 
a  factor  of  10. 


FRACTIONAL  BANDWIDTH  =  0.1 


Figure  3:  Threshold  coherence  value  from  (31)  v 

Gs/Gw  for  fractional  bandwidth  values  (7^)  (a) 
Gs/Gw  — >  00  are  superimposed  for  several  values 


FRACTIONAL  BANDWIDTH  =  0.5 


Let  us  examine  a  narrowband  signal  scenario  that  is  typical  in  aeroacoustics,  with  center  fre¬ 
quency  f0  =  lo0/(2tt)  =  50  Hz  and  bandwidth  A /  =  Auj/(2tt)  =  5  Hz,  so  the  fractional  bandwidth  is 
A f/fo  =  0.1.  From  Figure  3a,  coherence  |'ys|  =  0.8  requires  time-bandwidth  product  A f  -T  >  200, 
so  the  necessary  time  duration  T  =  40  sec  for  TDE  is  impractically  large  for  most  cases  with 
moving  sources. 

Increased  time-bandwidth  product  of  the  observed  signals  is  necessary  to  make  TDE  feasible  in 
environments  with  signal  coherence  loss.  As  discussed  with  respect  to  the  CRB,  only  the  observation 
time  is  controllable  in  passive  applications,  thus  leading  us  to  consider  source  motion  models  in 
Section  4  for  use  during  long  observation  intervals.  The  remainder  of  this  section  continues  to  focus 
on  non-moving  sources,  with  simulation  results  presented  in  Section  3.2  that  verify  the  CRB  and 
threshold  coherence  values  for  TDE,  and  a  discussion  in  Section  3.3  that  extends  the  H  =  2  sensor 
case  of  this  subsection  to  TDE  with  H  >  2  sensors. 

3.2  TDE  Simulation  Examples 

Consider  TDE  at  H  =  2  sensors  with  varying  signal  coherence  js.  Our  first  simulation  example 
involves  a  signal  with  reasonably  wide  bandwidth,  A /  =  30  Hz  centered  at  /o  =  100  Hz,  so  the 
fractional  bandwidth  A///o  =  0.3.  The  signal,  noise,  and  coherence  are  flat  over  the  frequency 
band,  with  SNR  Gs/Gw  =  100  (20  dB).  The  signals  and  noise  are  band-pass  Gaussian  random 
processes.  The  sampling  rate  in  the  simulation  is  Fs  =  101 2 3 4  samples/sec,  with  T  =  3  x  104  samples, 
so  the  time  interval  length  is  T  =  3  sec. 

Figure  4a  displays  the  simulated  RMS  error  on  TDE  for  0.2  <  js  <  1.0,  along  with  the 
corresponding  CRB  from  (26).  The  simulated  RMS  error  is  based  on  100  runs,  and  the  TDE  is 
estimated  from  the  location  of  the  maximum  of  the  cross-correlation  of  the  sensor  signals.  The 
threshold  coherence  for  this  case  is  0.41,  from  (31).  Note  in  Figure  4a  that  the  simulated  RMS 
error  on  TDE  diverges  sharply  from  the  CRB  very  near  to  the  threshold  coherence  value  of  0.41, 
illustrating  the  accuracy  of  the  analytical  threshold  coherence  in  (31). 

Next  we  consider  TDE  with  three  different  signals: 

1.  A  narrowband  signal  with  A /  =  2  Hz  centered  at  fo  =  40  Hz.  We  refer  to  this  as  “1 
harmonic”  (it  is  the  fundamental  frequency  of  the  signals  defined  next). 

2.  “2  harmonics”  at  40  and  80  Hz,  with  bandwidth  A /  =  2  Hz  at  each  harmonic. 

3.  “5  harmonics”  at  40,  80, 120, 160,  200  Hz,  with  bandwidth  A /  =  2  Hz  at  each  harmonic. 

The  signal,  noise,  and  coherence  are  flat  over  each  frequency  band,  with  SNR  Gs/Gw  =  100  (20 
dB),  and  the  signals  and  noise  are  band-pass  Gaussian  random  processes.  The  sampling  rate  in  the 
simulation  is  Fs  =  104  sanrples/sec,  with  T  =  2  x  104  samples,  so  the  time  interval  length  is  T  =  2 
sec. 

Figures  4b-d  display  the  simulated  RMS  error  on  TDE  (based  on  1,000  runs)  for  coherence 
values  0.7  <  <  1.0.  As  in  the  previous  example,  the  TDE  is  obtained  by  cross-correlation.  The 

threshold  coherence  is  defined  only  for  the  “1  harmonic”  signal,  and  the  threshold  coherence  value 
is  ~  1.  Figure  4b  illustrates  the  divergence  of  the  simulated  RMS  error  from  the  CRB,  except  at 
7  s  =  1- 

Figures  4c  and  d  display  the  results  for  the  signals  with  “2  harmonics”  and  “5  harmonics.”  The 
additional  harmonics  enable  accurate  TDE  for  lower  coherence  values,  but  we  cannot  use  (31)  to 
compute  the  analytical  threshold  coherence  for  the  harmonic  signals.  The  “approximate  threshold 
coherence”  values  indicated  in  Figures  4c  and  d  are  computed  as  follows.  For  the  K  harmonics, 
suppose  the  total  bandwidth  of  all  harmonics,  A /  =  I\  ■  2  Hz,  is  centered  at  the  fundamental 


frequency  /o  =  40  Hz.  The  approximate  threshold  coherence  values,  0.95  in  Figure  4c  and  0.50 
in  Figure  4d,  are  considerably  lower  than  the  actual  points  of  divergence  from  the  CRB.  Not 
surprisingly,  the  total  bandwidth  that  is  “spread”  across  the  harmonics  is  less  useful  for  overcoming 
signal  coherence  loss  than  an  equivalent  bandwidth  concentrated  at  the  fundamental  /q  =  40  Hz. 
Narrowband  and  harmonic  signals  are  generally  difficult  for  TDE  due  to  ambiguous  peaks  in  the 
cross-correlation  function. 


3.3  TDE  with  H  >  2  Sensors 


Weinstein  [14]  has  studied  TDE  with  a  network  of  H  >  2  sensors.  In  this  section,  we  extend  his  anal¬ 
ysis  to  partially  coherent  signals,  and  show  that  pairwise  TDE  is  essentially  optimum  for  cases  of 
interest  with  reasonable  signal  coherence  between  sensors.  By  pairwise  TDE  we  mean  that  one  sen¬ 
sor,  say  H ,  is  identified  as  the  reference,  and  only  the  H  —  1  time  differences  D\h ,  D2H ,  ■  ■  ■ ,  Dh~i,h 
are  estimated.  Under  the  conditions  described  below,  these  H  —  1  estimates  are  nearly  as  accurate 
for  source  localization  as  forming  all  pairs  of  TDEs  Dgh  for  all  g  <  h.  Weinstein’s  analysis  [14]  is 
valid  for  moving  as  well  as  non-nroving  sources. 

Assuming  equal  SNRC  at  all  sensors,  equal  coherence  between  all  sensor  pairs,  and  H-SNRC  3> 
1,  we  can  show  that  forming  all  TDE  pairs  Dgh  potentially  improves  the  source  localization  variance 
relative  to  pairwise  processing  by  the  factor 


(1  +  2-ijy) 


(32) 


Clearly  V  — ■>  1  as  7S  — ►  1,  and  V  <  (3iL)/[2(l  +  H)\  <  1.5  for  >  0.5.  Therefore  the  potential 
accuracy  gain  from  processing  all  sensor  pairs  is  negligible  when  the  coherence  exceeds  the  threshold 
values  that  are  typically  required  for  TDE. 

This  result  suggests  strategies  with  moderate  communication  bandwidth  that  potentially  achieve 
nearly  optimum  localization  performance.  The  reference  sensor,  H,  sends  its  raw  data  to  all  other 
sensors.  Those  sensors  h  =  1, . . . ,  H  —  1,  locally  estimate  the  time  differences  D\  h,  •  •  • ,  Dh~i,h, 
and  these  estimates  are  passed  to  the  fusion  center  for  localization  processing  with  the  bearing  es¬ 
timates  <(>1,...,  ■  A  modified  scheme  with  more  communication  bandwidth  but  more  centralized 

processing  is  for  all  H  sensors  to  communicate  their  data  to  the  fusion  center,  with  TDE  performed 
at  the  fusion  center. 


4  DATA  MODEL  FOR  A  MOVING  SOURCE 

Our  objective  in  this  paper  is  to  quantify  scenarios  in  which  jointly  processing  data  from  widely- 
spaced  sensors  has  the  potential  for  improved  source  localization  accuracy,  compared  with  incoher¬ 
ent  triangulation/tracking  of  bearing  estimates.  We  established  in  Section  2  that  the  potential  for 
improved  accuracy  depends  directly  on  TDE  between  the  sensors.  Then  we  showed  in  Section  3 
that  TDE  with  partially-coherent  signals  is  feasible  only  with  an  increased  time-bandwidth  product 
of  the  sensor  signals.  This  leads  to  a  constraint  on  the  minimum  observation  time,  T,  in  passive 
applications  where  the  signal  bandwidth  is  fixed.  If  the  source  is  moving,  then  approximating 
it  as  non-nroving  becomes  poorer  as  T  increases,  so  modeling  the  source  motion  becomes  more 
important. 

Approximate  bounds  are  developed  in  [15]  and  [16]  that  specify  conditions  of  validity  for  non- 
nroving  and  moving  source  models.  Let  us  consider  H  =  2  sensors  with  Doppler  values  a.2  >  «i  (see 


DELAY  (sec)  DELAY  (sec) 


RMS  TDE  (WIDEBAND):  f  =  100  Hz,  Af  =  30  Hz 


RMS  TDE:  1  HARMONIC,  f  =  40  Hz,  A  f  =  2  Hz 


(a) 


RMS  TDE:  2  HARMONICS,  f  =  40  Hz,  A  f  =  2  Hz 


(c) 


(b) 


RMS  TDE:  5  HARMONICS,  f  =  40  Hz,  A  f  =  2  Hz 


(d) 


Figure  4:  Comparison  of  simulated  RMS  error  for  TDE  with  CRBs  and  threshold  coherence  value, 
(a)  Wideband  signal  with  A /  =  30  Hz  centered  at  /o  =  100  Hz.  (b)-(d)  Narrowband  signal  with 
A /  =  2  Hz,  fundamental  frequency  /o  =  40  Hz  and  1,2,  and  5  harmonic  components,  respectively. 


(46)  for  the  definitions  of  ai,Q!2).  If  /max  (Hz)  is  the  maximum  signal  frequency  that  is  processed, 
then  TDE  estimation  accuracy  is  not  seriously  affected  by  ignoring  source  motion,  as  long  as  the 
time  interval  T  satisfies 

T< - -± - -.  (33) 


/max  (§*  -  l) 


Taking  typical  parameters  for  ground  vehicles  in  aeroacoustics,  let  us  consider  a  vehicle  moving 
at  5%  the  speed  of  sound  (15  m/sec),  with  radial  motion  that  is  in  opposite  directions  at  the  two 
sensors.  Then  012/oti  —  1  ~  0.1  and  (33)  becomes  T  <C  10//max.  For  /max  =  100  Hz,  the  requirement 
is  T  <C  0.1  sec,  which  according  to  the  analysis  in  Section  3.1  yields  insufficient  time-bandwidth 
product  for  partially  coherent  signals  that  are  typically  encountered.  Thus  motion  modeling  and 
Doppler  compensation  are  critical,  even  for  aeroacoustic  sources  that  move  more  slowly  than  in 
this  example. 

In  this  section,  we  extend  the  non-nroving  source  model  from  Section  2  using  first-order  motion 
models  (see  (34), (35), (49)).  The  first-order  motion  models  are  simple  and  accurate  over  larger 
time  intervals  T  compared  with  the  non-nroving  source  model.  However,  accurate  modeling  of 
more  complex  trajectories  over  longer  time  intervals  requires  higher-order  polynomial  models,  with 
added  complexity.  The  source  position  trajectory  is  modeled  as  a  straight  line  with  constant  velocity 
over  an  interval  of  length  T, 

xs{t)  =  xsfi  +  xs-  (t-t0),  t0<t<t0  +  T  (34) 

Vs(t)  =  Vsfi  +  Vs  ■  (35) 

so  xs,ys  are  the  velocity  components.  The  source  trajectory  parameter  vector  is 

©  =  [xs,o,xs,ysfi,ys]T,  (36) 

and  the  (time- varying)  propagation  time  from  the  source  to  the  sensors  on  array  h  follows  from  (1) 
and  (2): 


Th{t)  = 
Thn(t)  ~ 


dh(t) 


1  r 


C  L 


(xs(t)  -  xh )2  +  (ys(t)  -  yhf 


1/2 


dhit) 


dh{t) 


=  --  [(cos  <t>h(t))Axhn  +  (sin  (j)h(t))Ayhn\  ■ 

The  bearing  and  bearing  rate  are  related  to  the  source  motion  parameters  0  as 
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=  tan  1 

. ys(t )  -  yh. 

Vsfi  +  ys-(t~  to)  -  Vh 


4>h(t)  = 


ys  cos  4>h(t)  -  xs  sin  <ph(t) 
dhit) 


(37) 

(38) 

(39) 

(40) 


The  radial  velocity  of  the  source  with  respect  to  array  h  is 

vr,h(t )  =  cos  +  Vs  sin  (f>h(t).  (41) 

We  can  insert  (34)  and  (35)  into  (37)  to  obtain  the  following  approximation  for  the  propagation 
time  to  array  h: 


Th(t)  =  Th(t0) 


1  2  •  cos  (f>hjtp)  -Xs-jt-  to)  2  •  sin  < />h(t0)  •  ys  ■  (t  -  t0) 1 1/2 


dh(t0) 


dh (to) 


Th(to)  +  Vr'h ^  ■  (t  -  to), 


(42) 

(43) 


where  dh(to)  and  vr^{to)  are  the  source  distance  and  radial  velocity  at  the  start  of  time  interval 
t  =  to-  The  approximation  (43)  is  valid  as  long  as  the  total  motion  during  the  time  interval  T  is 
much  less  than  the  range,  i.e. ,  \2xsT\  <C  dh{to)  and  \2ysT\  <C  dh{to). 

Next  we  use  the  approximation  (43)  and  model  the  received  signal  at  the  reference  sensor  on 
array  h  as 

s h(t~Th(t))  =  sh  ^1  -  lr’h^to^  t,  -  Th(t0 )  +  Vr’h^to  (44) 

=  sh.  (ah  t  -  Th(t0 )  +  .  t0  <t  <t0  +  T,  (45) 

where 

ah  =  l-  =  1  -  i  [xs  cos  (/>h(to)  +  ys  sin  </>h(t0)]  (46) 

is  the  Doppler  compression  and 

Th{to)  =  dh ^  =  -  [(^,o  -  xh )2  +  (ys, o  -  yh)2}  1  (47) 

C  C  L 

is  the  propagation  delay  at  the  initial  time  t  =  to.  Without  loss  of  generality,  we  set  to  =  0,  so  the 
received  signal  at  sensor  n  on  array  h  is 

sh  (ah  t  -  rh( 0)  -  Thn(t))  ,  (48) 

which  is  the  extension  of  the  signal  component  of  (8)  to  the  moving  source  case.  Note  from  (2) 
that  t hn{t)  depends  on  the  source  location  only  through  the  time- varying  bearing  <j>h(t),  which  we 
approximate  with  a  first-order  model 

4>h{t)  «  <j)h(to)  +  <fr(to)  ■  (t  -  t0),  t0<t<to  +  T.  (49) 

For  a  single  array  h.  the  Doppler  compression  ah  and  time  delay  Th(to)  have  negligible  effect  on 
estimation  of  the  the  intra-array  delays  T^n(f),  since  ah  and  r^(fo)  are  identical  for  each  n  = 
1 ,  Nh-  Thus  each  array  can  be  processed  separately  to  estimate  the  bearings  (f>i(to), . .. ~ ,  <j>H(to ) 
and  bearing  rates  c/)i(to), . . .  ,<pH{to),  and  these  can  be  “triangulated”  via  (39)  and  (40)  to  estimate 
the  source  motion  parameters  ©  in  (36).  An  algorithm  [17]  for  estimating  (j>h(to)  and  <t>h(to)  is 
described  in  Section  5  and  demonstrated  with  measured  data  in  Section  6. 

Let  us  consider  the  signals  received  at  the  reference  sensors  at  each  array,  so  Thn{t )  =  0  in  (48): 

Sl  ^1  -  t  _  Tl^0)  ^1  -  Vr'H(to^  t  -  TH(to)  .  (50) 

Our  modeling  assumptions  imply  that  each  signal  Sh  ^1  —  ^  _  Th(f^  is  a  wide-sense  sta¬ 

tionary  Gaussian  random  process.  However,  for  two  arrays  g,h  with  unequal  Doppler  vr_g(to)  / 
Vr.h(to),  the  signals  at  arrays  g,h  are  not  jointly  wide-sense  stationary  [15,  18],  complicating  the 
analytical  description  and  the  CRB  performance  analysis.  The  jointly  nonstationary  sensor  signals 
generally  are  not  characterized  by  a  cross-spectral  density  matrix,  so  the  CRB  is  not  the  inverse  of 
a  FIM  of  the  form  (15).  An  approximate  CRB  analysis  for  TDE  with  jointly  nonstationary  signals 
as  in  (50)  is  given  in  [15].  The  CRB  analysis  is  rigorously  justified  in  [18]  and  extended  to  CRBs  on 
differential  Doppler.  A  clever  transformation  is  used  in  [18]  so  that  the  jointly  nonstationary  signals 
in  (50)  are  locally  modeled  by  a  CSD  of  the  form  (14),  and  it  is  shown  that  the  representation  is 
accurate  for  CRB  analysis. 


We  can  formulate  the  results  in  [18]  for  the  case  of  partially  coherent  signals4,  leading  to  the 
following  for  H  =  2  arrays,  assuming  large  tinre-bandwidth  product  (much  larger  than  the  coherence 
time  of  the  signals  and  noise).  We  define  the  TDE  D\2  =  T\  (to)  —  T2(fo)  and  the  differential  Doppler 
Avi2  =  vr,i(t0)  -  vr,2(t0). 

•  Estimation  of  TDE  and  differential  Doppler  are  decoupled,  so  the  CRB  on  D\2  is  given  by 
(25),  which  is  identical  to  the  non-nroving  source  case. 


•  The  threshold  coherence  analysis  for  TDE  in  (31)  and  Figure  3  extends  to  the  moving  source 
case.  In  the  best  case  that  Doppler  effects  are  perfectly  estimated  and  compensated,  the 
TDE  problem  that  remains  is  identical  to  the  non-nroving  source  case.  Doppler  estimation 
is  less  demanding  in  terms  of  tinre-bandwidth  product  compared  with  TDE.  Indeed,  Doppler 
estimation  is  possible  with  sinusoidal  signals  [18]  that  have  negligible  bandwidth. 


•  The  CRB  on  differential  Doppler  [18],  modified  for  partially-coherent  signals  and  assuming 
SNRC!i(w)  =  SNRc,2(o;)  =  SNRc(w),  is 


CRB(Aui2) 


24tt  f  c\2  [o  fu»  io2  SNRc(u;)2  ^  ' 

~T  yr)  2  Jo  1  +  2  •  SNRc(ca)  duJ 


(51) 


Note  that  (51)  is  a  scalar  multiple  of  the  CRB  on  TDE  in  (25).  However,  the  CRB  on 
differential  Doppler  may  be  achievable  in  scenarios  in  which  the  tinre-bandwidth  product  is 
insufficient  for  TDE. 


Interestingly,  differential  Doppler  provides  sufficient  information  for  source  localization,  even 
without  TDE,  as  long  as  five  or  more  sensors  are  available  [18].  Thus  the  source  motion  may 
be  exploited  in  scenarios  where  TDE  is  not  feasible,  such  as  narrowband  signals  [18]. 


•  We  discussed  TDE  with  H  >  2  sensors  in  Section  3.3,  concluding  that  pairwise  processing  of 
TDEs  D ih,  . . . ,  Dh-i,h  with  a  reference  sensor  H  is  nearly  optimum  for  scenarios  of  interest. 
A  similar  result  holds  for  differential  Doppler  estimation  [18],  where  pairwise  estimation  of 
Avm,  ■  ■  ■ ,  Avh~i,h  is  nearly  as  accurate  as  estimation  of  all  pairs  A vgh,  as  long  as  H  • 
SNRc(w)  >  1. 


5  AN  ALGORITHM 

The  parameters  that  can  be  directly  estimated  from  the  sensor  data  are  the  bearings  (f>i(to), . . ., 
0#(to)>  bearing  rates  <j>i(to), . . .  ,<j>H{to),  pairwise  time  differences  Dm  =  t\  (to)  —  T#(io),  ..., 
DH-i,h  =  TH-i(to)  -  TH(t0),  and  differential  Doppler  A v1H  =  Wi^o)  -  vr,H{t0), . . . ,  AvH-i,h  = 
vr,H-i(to)  —  v-r.HiU))-  Equations  (39),  (40),  (46),  (47)  define  the  nonlinear  relations  that  “triangu¬ 
late”  these  parameters  and  relate  them  to  the  source  motion  parameters  ©  =  [xSio>  xS:  ys^,  ys]T ■ 

A  distributed  processing  algorithm  is  outlined  below,  and  parts  of  the  algorithm  are  illustrated 
with  measured  aeroacoustic  data  in  the  next  section. 

1.  Use  the  local  polynomial  approximation  (LPA)  beamformer  [17]  at  each  array  to  estimate 
the  bearings  and  bearing  rates.  The  LPA  beamformer  in  [17]  is  formulated  for  narrowband 
processing,  and  it  is  a  generalization  of  the  classical  beamformer  to  moving  sources.  We 
extend  it  in  a  straightforward  way  to  wideband  signals  by  incoherently  averaging  the  LPA 
beampatterns  at  different  frequencies. 

4The  signal  coherence  between  the  signals  at  arrays  g  and  h  in  (50)  is  defined  assuming  perfect  compensation  of 
the  Doppler  compression  ag,cth,  thus  yielding  the  definition  in  (7). 


2.  Solve  (39),  (40)  to  obtain  initial  estimates  of  the  source  motion  parameters  0.  These  estimates 
correspond  to  incoherent  triangulation  of  the  bearings  and  bearing  rates  from  individual 
arrays. 

3.  Estimate  the  Doppler  compression  factors  a\, . . .  ,an,  compensate  for  Doppler,  and  test 
whether  the  signals  at  distinct  arrays  have  sufficient  coherence,  fractional  bandwidth,  and 
time-bandwidth  product  to  enable  TDE  between  arrays  (see  Section  3.1  for  the  conditions). 

4.  If  the  conditions  are  not  met,  then  incoherent  triangulation  of  the  bearings  and  bearing  rates 
is  nearly  optimum,  and  further  joint  processing  is  not  informative. 

5.  If  the  conditions  are  met,  then  identify  a  reference  array  H  (the  array  with  maximum  SNR) 
and  estimate  the  time  differences  Dm,---,  Dh~i,h  and  differential  Dopplers  Avm,  •  ••, 
A vh-i,h-  Many  estimation  methods  have  been  studied  for  this  problem,  e.g.,  [19]- [27] .  The 
maximum  likelihood  solution  involves  wideband  ambiguity  function  search  over  Doppler  and 
TDE  [15],  while  the  deskewed  short-time  correlator  [22]  is  a  computationally  simpler  approx¬ 
imation. 

6.  A  suboptimum  procedure  is  to  avoid  the  joint  Doppler  and  TDE  estimation  in  the  preceding 
step,  and  instead  use  the  initial  Doppler  estimates  from  steps  1  and  2  and  perform  TDE  after 
approximate  Doppler  compensation.  With  this  approach,  triangulation  of  the  TDEs  via  (47) 
will  improve  the  estimates  of  xSjo  and  ys> o  only  (and  not  the  source  velocity  xs,ys). 

7.  If  multiple  sources  are  present,  then  the  LPA  beamformer  in  step  1  may  be  used  to  separate 
the  source  signals  at  each  array  prior  to  Doppler/TDE  estimation.  We  note  that  wideband 
beamforming  algorithms  for  non-nroving  source  models  are  presented  in  [28,  29]  and  [30]- [32]. 

The  LPA  beamformer  in  steps  1  and  7  is  illustrated  in  the  next  section  for  a  two-source  scenario 
based  on  measured  aeroacoustic  data.  Examples  of  TDE  with  Doppler  compensation  (step  5)  are 
presented  in  [33]. 

6  EXAMPLE  WITH  TWO  MOVING  SOURCES 

We  present  an  application  of  the  local  polynomial  approximation  (LPA)  beamformer  [17]  to  mea¬ 
sured  aeroacoustic  data  with  two  ground  vehicles:  Ml  and  M3  tanks.  The  vehicle  trajectories  over 
a  10  second  segment  are  shown  in  Figure  5a.  Two  arrays,  labeled  8B  and  8C,  are  separated  by 
about  23  nr.  Each  array  is  circular  with  N  =  7  sensors,  4-ft  radius,  and  six  sensors  equally  spaced 
around  the  perimeter  with  one  sensor  in  the  center.  We  present  results  based  on  processing  the 
data  at  array  8B,  for  which  the  bearing  and  bearing  rate  of  the  sources  are  shown  in  Figure  5b. 
The  bearing  of  the  M3  varies  by  more  than  30°  over  the  10  second  time  interval.  The  range  of  the 
M3  is  closer  at  about  100  nr,  with  the  Ml  range  at  approximately  200  nr. 

The  data  from  array  8B  is  processed  over  a  wide  bandwidth  from  30  to  150  Hz.  The  beanrpat- 
tern  for  a  classical  beamformer  based  on  a  non-nroving  source  model  is  shown  in  Figure  5c.  The 
beanrpattern  is  the  incoherent  sum  of  narrowband  beanrpatterns  over  the  30  to  150  Hz  frequency 
band.  The  peak  of  the  beanrpattern  is  located  at  approximately  the  mean  bearing  of  the  stronger 
source  (M3)  over  the  10  second  interval. 

The  beanrpattern  of  the  LPA  beamformer  is  shown  in  Figure  6a.  The  LPA  beanrpattern  exploits 
a  first-order  model  for  time- varying  bearing,  as  in  (49).  The  LPA  beanrpattern  is  two-dinrensional, 
with  axes  of  initial  bearing  4>(to)  and  and  bearing  rate  (j>(to).  The  LPA  beamformer  in  [17]  is 


formulated  for  narrowband  processing,  so  we  extend  to  the  wideband  case  by  incoherently  adding 
the  narrowband  LPA  beampatterns  over  the  frequency  range.  The  peak  of  the  beampattern  in 
Figure  6a  is  close  to  the  true  values  of  (j>(to),  4>{to)  for  the  M3  that  are  shown  in  Figure  5b.  The 
location  of  the  weaker  Ml  source  is  not  evident  in  the  LPA  beampattern  in  Figure  6a,  so  we  subtract 
an  estimate  of  the  stronger  source  from  the  data.  The  subtraction  is  performed  based  on  the 
bearing  and  bearing  rate  estimates  from  Figure  6a.  The  subtraction  is  coherent  over  the  processing 
bandwidth  and  includes  the  time- varying  bearing.  The  LPA  beampattern  after  subtraction  is  shown 
in  Figure  6b,  which  indicates  the  bearing  and  bearing  rate  of  the  weaker  Ml  source.  This  example 
illustrates  the  gain  in  resolution  that  is  achieved  by  exploiting  the  source  motion  in  a  beamformer. 

7  CONCLUDING  REMARKS 

The  potential  gain  in  source  localization  accuracy  when  data  from  distributed  arrays  is  processed 
jointly  and  coherently  is  quantified  by  the  CRBs  presented  in  this  paper.  The  amount  of  improve¬ 
ment  and  the  feasibility  of  achieving  the  improvement  depend  critically  on  the  scenario,  which  is 
characterized  by  the  coherence  between  source  signals  arriving  at  distributed  sensors,  the  signal 
bandwidth  and  spectrum  shape  (wideband  vs.  harmonic),  the  observation  time  for  coherent  pro¬ 
cessing,  the  noise  level,  the  source  motion  parameters  (velocity,  complexity  of  maneuvers),  and  the 
number  of  sources.  In  feasible  scenarios  in  which  the  time-bandwidth  product  is  large  enough  to 
enable  TDE,  we  presented  an  algorithm  that  requires  moderate  communication  bandwidth  between 
sensors.  The  processing  involves  estimation  of  bearing  and  bearing  rate  at  individual  arrays,  and 
estimation  of  time  delay  and  differential  Doppler  between  pairs  of  arrays. 


References 

[1]  R.R.  Tenney  and  J.R.  Delaney,  “A  distributed  aeroacoustic  tracking  algorithm,”  Proc.  American  Con¬ 
trol  Conf.,  pp.  1440-1450,  June  1984. 

[2]  Y.  Bar-Shalom  and  X.-R.  Li,  Multitarget-Multisensor  Tracking:  Principles  and  Techniques,  YBS,  1995. 

[3]  A.  Farina,  “Target  tracking  with  bearings-only  measurements,”  Signal  Processing,  vol.  78,  pp.  61-78, 
1999. 

[4]  B.  Ristic,  S.  Arulampalam,  C.  Musso,  “The  influence  of  communication  bandwidth  on  target  tracking 
with  angle  only  measurements  from  two  platforms,”  Signal  Processing,  vol.  81,  pp.  1801-1811,  2001. 

[5]  L.M.  Kaplan,  P.  Molnar,  Q.  Le,  “Bearings-only  target  localization  for  an  acoustical  unattended  ground 
sensor  network,”  Proc.  SPIE  AeroSense,  Orlando,  Florida,  April  2001. 

[6]  R.J.  Kozick  and  B.M.  Sadler,  “Distributed  Sensor  Array  Processing  of  Wideband  Acoustic  Signals,” 
1999  Meeting  of  the  IRIS  Specialty  Group  on  Battlefield  Acoustics  and  Seismics,  Laurel,  MD,  September 
13-15,  1999. 

[7]  R.J.  Kozick  and  B.M.  Sadler,  “Algorithms  for  Localization  and  Tracking  of  Acoustic  Sources  with 
Widely  Separated  Sensors,”  Proc.  2000  Meeting  of  the  MSS  Specialty  Group  on  Battlefield  Acoustics 
and  Seismics,  Laurel,  MD,  October  18-20,  2000. 

[8]  B.  Friedlander,  “On  the  Cramer-Rao  Bound  for  Time  Delay  and  Doppler  Estimation,”  IEEE  Trans, 
on  Info.  Theory,  vol.  IT-30,  no.  3,  pp.  575-580,  May  1984. 


NORTH  (m) 


VEHICLE  PATHS  AND  ARRAY  LOCATIONS:  710  to  720  SEC 


120  - 

100  - 

80  - 

60  - 


50 

EAST (m) 


(a) 


■•+•  Ml 
x  M3 

□  ARRAY  8 B 
O  ARRAY  8 C 


LU  3, 

!< 

IT 

(5  2 

Z 

tr 

ft  1 


Ml 

M3 


714  715  716 

TIME  (sec) 

BEARING  RATE  B 


. x’ ' 

1  O  Ml 

|  -  x  M3 

O . 

"O . O . o 

- 

. o 

712  713  714  715  716  717 

TIME  (sec) 


(b) 


719  720 


CLASSICAL  BEAMFORMER  B 


Figure  5:  (a)  Array  locations  and  trajectory  of  Ml  and  M3  vehicles  over  a  10  second  time  interval, 
(b)  Bearing  and  bearing  rate  of  sources  with  respect  to  array  8  B.  (c)  Beampattern  of  classical 
beamformer  based  on  non-nroving  source  model. 
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Figure  6:  (a)  Beampattern  of  LPA  beamformer,  showing  the  location  of  the  stronger  source  (M3), 
(b)  Beampattern  of  LPA  beamformer  after  subtracting  an  estimate  of  the  strong  source  signal  from 
the  data,  showing  the  location  of  the  weaker  source  (Ml). 
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